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Abstract 

We analytically compute a localization criterion in double scattering ap- 
proximation for a set of dielectric spheres or perfectly conducting disks uni- 
formly distributed in a spatial volume which can be either spherical or layered. 
For every disordered medium, we numerically investigate a localization crite- 
rion, and examine the influence of the system parameters on the wavelength 
localization domains. 
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I. INTRODUCTION 



Several works of Anderson on localization |IJ have mostly been related to electrons trans- 
port in solids with random impureties. Wolfle has proposed a diagrammatical treatment 
for the electronic localization in a bidimensional disordered medium. This formalism was 
extended to the electromagnetic waves in a volume of disordered medium [fj, Q], where 
the impureties are modeled by dielectric spherical scatterers, and the role played by the 
backscattering mechanism in the energy localization was established. The physical basis of 
the localization is a consequence of the vanishing of diffusion coefficient of the electromag- 
netic wave energy, this condition was then called a localization criterion. 

The purpose of this paper is to analytically determine several conditions on the param- 
eters of some disordered systems (namely, the number of scatterers, their size, permittivity, 
geometry, the width of the medium) satisfying a localization criterion of the electromagnetic 
energy, and to show their influence on the wavelength localization domains ||. We give 
the analytic expressions of the localization criterion for a spherical or layered medium made 
of dielectric spheres or perfectly conducting disks. We numerically determine the sets of 
parameters offering ranges in wavelength where the localization criterion is satisfied. More 
precisely, we compute the evolution of the range when we modify every parameter, and we 
show that the location and the bandwitdh of these ranges may be controlled by varying 
some characteristic parameters of the disordered system. 

This paper is organized as follows. After the introduction, we recall in Section 2 the 
mathematical model used to describe the behaviour of the electromagnetic wave in a disor- 
dered medium, and the expression of a localization criterion. Then in Section 3, we derive 
analytically a localization criterion for several disordered media. In Section 4, we numeri- 
cally compute the localization criterion for these media in order to determine localization 
domains, and we examine the influence of the system parameters on these localization do- 
mains. Computations with a multidipolar formulation are used for comparison with the 
theoretical model. Finally, Section 5 is devoted to the conclusions. 



1 



II. THE MATHEMATICAL MODEL 



We consider N identical scatterers of relative permittivity e s (u>) distributed in a spatial 
volume V. The location of scatterers is given by the vector R p ,p G [1, N]. The background 
medium has a relative permittivity eo- An incident electromagnetic plane wave of pulsation 
uj, and wavevector k = — y/eo, interacts with the set of scatterers (figure [T]). 
The electric field E{y*, uj) at location ~r* satisfies the Helmholtz equation 

Ai^T*, uj) + k 2 e(-r>, w).T2(T>, uj) = 0, (1) 

where k = — is the wavevector in the vacuum, e(T*,uj) is the second rank tensor associated 

c 

with the relative permittivity of the medium, where the components are given by 

N 

eij(T>,uj) = eodij + £ (e sij (uj) - eo^Xj^r*), (2) 
p=i 

is the characteristic function of the p th scatterer 



1 if T* inside the scatterer, 
\/t;;(^)= { (3) 
outside. 

We assume the scatterers homogeneous and isotropic, thus 

ey (!*£<;)= ^e + ^(e s M -e )x^(T > )^ ii , (4) 

so we can treat the relative permittivity of the medium as a function of and uj. 
Then, we write the Helmholtz equation in a perturbative form 

AE*(T>,uj) + [kl-5V{J>)YE!{J>,uj) =0, (5) 

where 5V(~r*) = k 2 (eo—e s ) J2iLi Xj^(^) is the potential function describing the perturbation 
of the incident wavevector. 
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A. The Green's formalism 



In order to solve the Helmholtz Eq. (5), we first introduce the dyadic Green function 
gij(~r*, r, to) || which determines the wave at location ~r* generated by a pointlike cur- 
rent source located in r . Its components satisfy the Helmholtz equation 



Agij(~r > , r ,u) + [kl - 5V(T > )J^-(T > , r' ,u) = 5 ij -5(T > - /). (6) 

The interest of introducing a dyadic Green function is to provide a generic formulation for 

the electromagnetic field behaviour. 

This equation can be written in an integral form 

- — 3 f - — 

ft i (T > , T",u) =5-?-(T > -T >/ ,o;) + ^ / 9vSj^ ~ rt,v)5V(rt)g kj (ri, T >/ , u)dft , 

k=\ J 

(7) 

where <^(~r* — r,u) is the dyadic Green function of the electromagnetic wave in the non- 
perturbated medium of relative permittivity e , which satisfies the equation 

- = (Sa + ftSiV ( r - r"), (8) 

and g°(~f ¥ — T* 7 ) is the scalar Green function of the wave 

g\T> = -- * (9) 

An \ r — r'\ 

We introduce the Green operators Q and Q° 

G° : ~&(r) ^ J g =0 (r -r')~&(r')dr\ (10) 
: ~0?(r) -»• / g(r,r'Y${r')dr'. (11) 



The dyadic Green functions are defined by the kernels of the operators. The operator 5V 
in Eq. (7) is defined as the multiplicative operator, and the direct product between two 
successive operators is written point. 
Using this formalism, we can rewrite the integral equation like 

g = g° + g°.5V.g. (12) 



B. First moment of Q - Self energy operator 

1. Model of disordered medium 

The calculation of the kernel of Q using a perturbative expansion of (12) requires the 
knowledge of the 3iV space coordinates relative to the location of scatterers. Due to their 
complexity, we introduce a statistical approach of the problem. We call pn(Ri---Rn) the 
probability density of location of the N scatterers in the medium M , and we introduce the 
mathematical expection value or mean IE defined by 

: ?(T>) -> f 9ij (-T - r >, ,a;)p iV (^...^)nf =1 ^ p V(^ / )^ , > (13) 
where we average over the whole possible configurations of the disordered medium. 



2. First moment of Q 

From Eq. (12), (13) we obtain 

e(s) =g° + g°~E(5v.g). (14) 

By introducing the self-energy operator E J|, defined by 

E = ]E(SV.g^ [e((/)] -1 , (15) 
we can rewrite (14) like a Dyson equation 

e(^) = g° + g°TM(g). (16) 

The definition (15) of E is not in a convenient form to obtain a solution, because E^C/J is the 
unknown of the equation. So, in order to calculate E, we prefer to link it with the multiple 
scattering formalism ||. If T is the scattering operator of the medium, T is defined by 

g = g° + g°.T.g°, (17) 
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and the expectation value E(T) of the scattering operator is given by 

e(T) = g° + g°.iE(T).g . (is) 

We have the following relation between S and E(T), 



E = E(T) 



n -i 



1 + £°.E(T) 



E(T) - E(T).C/°.E(T) + (19) 



where T can be expressed using a multiple scattering expansion in terms of the scattering 
operator by a single scatterer ti at location R-i, 

N N N N N N 

r = $> + E*i0°- E *; + E*iS°- E E ** + ■■■ (20) 

i=l i=l j=\ : jjti i=\ j=ijjLi k=l,ky£j 

The interest of the last expression is to express S by the mean of a scattering theory instead 
of calculating it from a perturbative expansion. 

3. Macroscopic homogenization of the disordered medium - Scalar formulation 

We use the hypothesis of macroscopic homogenization of the disordered medium in order 
to calculate the kernel of E^C?), i.e., the dyadic Green function of the averaged wave. We 
suppose that after averaging, the medium behaves like a homogeneous one, which implies 
an invariance under translation. In consequence, we can first write the Dyson equation (16) 
like a convolution equation 

E(S) =£° + £°*X*E(V). (21) 

The other consequence of the macroscopic homogenization hypothesis is the possibility to 
write a scalar form of the previous equation. Indeed, the averaged medium being homoge- 
neous, the kernel of E^C?) may be written as 

E(^(T > - t')).. = (Sij + did^EfaCT - T")), (22) 

where E^("r > — "r*')) is the first moment of the scalar Green function of the electromagnetic 
wave. The same consideration holds for the kernel of the self-energy operator. Extracting 



the scalar part of Eq. (21), we obtain the scalar expression of the Green function in Fourier 
space 



E ( G(T ^ K^-W,,) - (23) 



where we notice that the invariance by translation of the averaged medium implies that the 
first moment of the Green function and the self energy function in Fourier space are only 
k dependent. E(G) is the scalar Green function of the homogenized medium. 

The argument of Sheng [|| is to treat the averaged medium like an effective medium, 
where the characteristic length of the averaged inhomogenities is small with respect to 
the wavelength, and therefore, they are not resolved by the wave, which explains that the 
averaged medium may be considered like a homogeneous one. In that case, the self-energy 
function £ describing the local microstructures of the medium, is weakly dependent on the 
variable k , so it is treated like a pulsation dependent function E(o»). 

If this condition is satisfied, the first moment of the Green function physically represents 
a propagative wave of effective wavevector k\ = k\ + S(u;), from which we deduce the 
wavevector of the wave 



k e = Re (fc e ) 



I + 4/I 
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(24) 



where T = Re^S(u;)^, 7 = Im^£(cG>)), and the scattering length I [[| which is the inverse of 
twice the damping rate (3 of the modulus of the averaged wave 

' = £ = v (25) 

The decreasing factor in the expression of the first moment of the Green function shows 
the fact that during the path in the averaged medium, the wave loses its phase coherence 
over a characteristic length given by Z, as a consequence of the successive scatterings. Over 
a distance of several scattering lengths I, the first moment of the Green function becomes 
negligible. In order to describe the behaviour of the electromagnetic wave for large distances, 
we need to investigate the second moment of the dyadic Green function. 
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C. The second moment of Q - Localization of the energy 



The averaged energy operator V is given by the averaged tensorial product of the two 



Green operators Q and Q* , where their kernels are respectively the dyadic Green function 

5to 

of the electromagnetic wave at pulsation oj + — uj + — and the conjugate complex dyadic 
Green function at pulsation to^ = u — 



e(p) = E(p®g* 



(26) 



where to is the central pulsation of the two dyadic Green functions, and Suj is their pulsation 
difference. We notice that 8u is the conjugated variable of the time traject t of the electro- 
magnetic wave in the medium. The behaviour for t — ► +00 of the averaged energy will be 
given by the behaviour of the kernel of TE(V^J when 5 to —>■ 0. 

The outer product ® applied on the components of the dyadic Green functions is given by 



(27) 



Using the macroscopic homogenization hypothesis, we can write the scalar form of the 
kernel of E(p), which is defined as the function e(p(T*, r' , to, Suj)\ and is related to the 
average energy density function L by its spatial Fourier transform 



e(p(^,?",^,^)) = — / L{~t lUj^uy^-^-^d-t . 

^ ' 2,71 J 



(28) 



Let us remark that the q variable in the function L is conjugated with the wave traject 
IT* — ~~r*'\. The behaviour of the average energy at large distance will be given by the 
behaviour of L when q — > 0. 

Following the frameworks of Arya || and Sheng ||, at large distances (q — > 0) and large 
time (Sou — > 0), L is governed by a diffusion-like equation 



—iSuL(~q y , to, 5uj) = q 2 D(~q* ', uj, 5u)L(lf ', to, 5to), 



where 



71 



3kll 2 Jq=o 



Q 2 



iSu 



-dQ 



(29) 



(30) 



is the general energy diffusion coefficient of the wave, and 




ujI 
3k e 



(31) 



is the Boltzmann diffusion coefficient. The corrective term involved in the general expression 
of the diffusion coefficient is a consequence of the backscattering effect of the electromagnetic 
wave in the disordered medium. 

The energy of the electromagnetic wave is said to be localized when the general diffusion co- 
efficient (30) vanishes, i.e. the corrective term generated by the backscattering compensates 
the Boltzmann diffusion coefficient, which gives us the condition 



k e and / are given by the relations (24) and (25). This condition is called the localization 
criterion of the electromagnetic wave energy. 



We propose to derive an analytic expression of the localization criterion for an electromag- 

27T 

netic wave of wavelength Ao = — interacting with the two following sets of media. 

k 

- A set of N identical dielectric spheres of radius a and relative permittivity e s , uniformly 
distributed in a spherical volume V , or in a layered volume of thickness e z (figure 0). 

- A set of N identical perfectly conducting disks of radius a and relative permittivity e s , 
uniformly distributed in a spherical volume V , or in a layered volume of thickness e z (figure 



The definitions of the averaged wavevector k e and the scattering length I are given by the 
relations (24) and (25). To express the k e .l product, we have to calculate the self-energy 
function for every disordered medium. 




(32) 



III. CALCULATION OF THE LOCALIZATION CRITERION 



FOR SEVERAL MEDIA 



1 



S 



A. Macroscopic homogenization conditions 



First, we need to verify some conditions on the media in order to consider them as macroscop- 
ically homogeneous, which means that the self-energy function must be only k dependent. 

The self-energy function in Fourier space for a medium of N identical scatterers of scat- 
tering function t( k , k ') is given by 

E("jfe\ t') = yt(t, t') J e^-^'^dlti, (33) 

where V is the volume of the medium. The integral in (33) taken over a finite volume can 
be approximated as 



J^-^d% = 60t-^ if< 



radius R of the medium >> — 

A 4 (34) 

thickness e z » — 
2 

which leads to the condition of macroscopic homogenization, because the delta distribution 
argument k ' — k implies that the self-energy function is only k dependent. 



B. Calculation of k e and I 

Using expansions (19) and (20), the first and the second scattering order of the self energy 
operator are respectively given by the following expressions 

N 

E = E(5>), (35) 
i=i 

N N 

E = E(J>) -E(5>-0°.*i)- (36) 
i=i i=i 

By assuming the macroscopic homogenization hypothesis, we can derive the scalar expres- 
sions of the kernels. For the second scattering order, we suppose that every scatterer receives 
a wave being locally plane, in that case, the scalar part of the Green propagator simplifies to 
<7°( k ,u) — —i-K5{k 2 — Since the iV scatterers are uniformly distributed in the volume 
V, we obtain 

N 

£(c<-0 = 'yt( U} ) sm §l e scattering order, (37) 



L(u;) = — t(a>) + 7r — -t ltd) second scattering order, 
V V 



(38) 



where 



V = I 



|7ri? 3 for a spherical medium of radius R 
e x .e y .e z for a layered volume of thickness e z 
and transverse widths e T and e„ 



(39) 



R or e x , and e z must satisfy the conditions (34). The scattering functions are given by 



t{u) 



ffE^(2p+l)(a P + M 



for a sphere of radius a and permittivity e s \ 

- kid^i^f) 2 + 1 (^) 4 + *f?(¥ 



(40) 



for a little disk of radius a (d =2a) |5|. 
From (38), (39), (40), we obtain the expression of the averaged wavevector k e and the 
scattering length / at the second scattering order. 



- For iV spheres uniformly distibuted in a volume V 



k e — p= 



1 + 



ft Th / \ 2\ 



P =i 



P =i 



L / | Eg°=i( 2 P + IJjMfh + M + t (MEgi(2p + l)(a p + & P 
^ A; 2 + f E~i( 2 P + l)Re(a p + b p ) + f (lm(E+^(2p + l)(a p + b p )) 



and 



(41) 



f E p °°=i(2p + 1)Mop + M + x( Im (E P + = °?(2p = l)(op + 6p)) 



2 ' 



(42) 



For N disks uniformly distributed in a volume V 
1 



V2 
1 + 



ifcg - 2nv / 2vrA; 2 l dRe( J D§) + 8Tm 2 d?k$(Re(D%) 2 



1 + 



\ y k 2 Q - 2nV2^k 2 Q dRe(D° ) + 87tn 2 d 2 k%((Re(D° 



(43) 
(44) 
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< 45 > 



and 

/ = * 

-2nV27rk 2 dIm(D$) + 87m 2 d 2 k$(lm(D 
where D° = f (^f ) 2 + ^§( h ff + ^(^fY, and n = N/V is the density. 

IV. NUMERICAL RESULTS 

A. Characterization of the localization domains 

We are now in position to search for localization domains, i.e, ranges of wavelength 

R 

satisfying the localization criterion k e .l < \ —. In the following we have computed the k e .l 

V 7T 

product as a function of the incident wavelength A, and labeled the localization curve as 
k e .l(X). Then we will look for the influence of the system parameters on the localization, 
for instance: the number N, the width a of the scatterers and for the case of spheres, their 
permittivity e s , the volume V of the medium. In our simulations, we have supposed that 
the relative permittivity e s is weakly dependent on the pulsation, so we have treated it as a 
scalar. However, it is possible to give a description by a frequency law dependent upon the 
nature of the scatterer. The calculations were performed using Mathematica . 

The first medium we study is an array of iV dielectric spheres of radius a, permittivity 
e s , uniformly distributed in a volume V of radius R. We numerically found a localization 
domain for the following set of parameters (figure |j) 
N = 1 million 
a = 0.01// 
e s =16 
R = 2fx 

A E [0.0805/x, 0.083//]. 

We satisfy the criteria of macroscopic homogenization of the medium because in that case 

R R 1 

the ratio — is about 25, and so the condition — >> - is valid. Let us notice that we are 
A A 4 
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working in the range k.a < 1, and only the first term of the expansion (40) was used in the 
calcultations of the scattering function for the spheres. 

Now, we will investigate the effect of the parameters on the localization domain. 

We show in figure [5], the localization curves for several values of the spheres radius. The 
dotted curve is taken for reference (a = O.Olyu). When the value of a decreases (0.095/x) 
(plain curve), the minimum of the curve occurs for lower wavelength, the corresponding 
value of the k e .l product increases and the localization domain narrows. Oppositely, when 
the radius is increased (0.0105/x) (dashed curve), the minimum of the localization curve 
occurs for higher wavelength, the corresponding value of the k e .l product decreases and the 
localization domain enlarges. 

In figure |], we start from the reference configuration, and we modify the relative permit- 
tivity of spheres. The curves are shown for the values: e = 15 (plain curve), e = 16 (dotted 
curve) and e = 17 (dashed curve). When we slightly increase the relative permittivity of the 
spheres, the localization domain enlarges with a translation to higher wavelengths. 

Next, we have fixed the real part of the permittivity (Re(e) = 16), and we simulate a 
absorption on the surface of spheres by adding an imaginary part to the relative permittivity, 
so we obtain a dissipative dielectric. The figure [7] shows the localization curves for the 
respective values (plain curve), 0.4 (dotted curve), 0.8 (dashed curve) and 1.2 (long dashed 
curve) of the imaginary part of the relative permittivity. When the value of the imaginary 
part increases, the value of the minimum of the localization curve also increases and the 
localization domain narrows and finally disappear. By adding an imaginary part to the 
relative permittivity, the scatterers absorb a part of the electromagnetic energy and the 
interactions during successive scatterings are reduced, which explains the disappearance of 
localization. 

N 

An other interesting parameter is the density n = —. Indeed, when we increase the 
density of scatterers, we also increase the backscattering contributions from the successive 
interactions of the wave with scatterers, as a consequence the localization is more easily 
attained. 
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We verify the influence upon the number N of scatterers in a fixed volume medium. 
The Figure |8| shows the localization curves corresponding to N — 1 million (plain curve), 
N —1.5 millions (dotted curve) and iV =2 millions (dash curve) spheres. When the number 
of scatterers increases, we observe a deepening of the localization curve minimum and a 
broadening of the localization domain. The inverse phenomenon was obtained when we fix 
the number of spheres and at the same time increase the radius of the medium. Figure |9| 
shows the localization curve for the radius R values 2, 2.1, 2.2, 2.3, and 2.4 /i respectively. 
In that case, the minimum of the localization curve increases with R and the localization 
domain disappears. The localization gradually vanishes because the medium is diluted, so 
the interactions between the scatterers are reduced. In both cases, the N or R variations 
show that the location of the minimums are weakly dependent, it means that one can choose 
a configuration giving a central vawelength localization, and then modify the width of the 
localization domain around this value by the mean of medium density. 

In the next part, we only focuse on the density effect. 

In the following example, we consider N dielectric spheres uniformly distributed in a 
rectangular volume [— ^f, y] x [— y, ^f] x [— y, y], where e x and e y are large in front of A 
in order to approximate the rectangular medium by a single layer of thickness e z . 



We numerically found a localization domain for the set of parameters (figure [H|) 

N= 10 5 
a = 0.01// 

thickness e z = 1.5 /i 
transverse lengths e x =e y =1.5 fi 
A G [0.0805^,0.083//]. 



The figure |TT] presents the localization curves when the number N of spheres takes 
respectively the values: 10 5 (plain curve), 1.210 5 (dotted curve) and 1.510 5 (dashed curve), 
the other parameters are kept fixed. We observe that the minimum of the curves diminishes 
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with the increase of N and the localization domain enlarges. Then, if we fix the number of 
spheres iV = 10 5 , and modify the thickness e z of the layer with respect to the macroscopic 



homogenization condition. The figure [12] shows the curves obtained for e z = 1/x (plain 
curve), e z = 1.5 /j, (dotted curve) and e z = 1.8// (dashed curve). The minimum of the curves 
increases with the thickness e z of the layer, the localization domain narrows and finally 
disappears. The observed effects are the same as in the case of a spherical volume. 

In a second series of tests, we replace the spheres by perfectly conducting disks, oriented 
following a plane perpendicular to the wavevector (figure |3[). 

For the two respective media (spherical and layered), we numerically find a localization 
domain for the set of parameters in a spherical volume (figure |l"3| ) , 
N = 1 billion 
a = 0.05/i 

volume radius R = 12.5// 

A E [0.45/1, 0.64/i] 

and for a single layer (figure ]14|) 

iV = 150millions 

a = 0.05 fi 

thickness e z = 10 \i 

transverse lengths e x =e y =10 \i 

A G [0.45//, 0.74//]. 



Figures [15] and [16] show the localization curves for the spherical volume case when we vary 
respectively the number of disks and the radius of the volume. In figure [15], the volume radius 
is fixed, R = 15//, and the number iV takes the values 1 billion (plain curve), 1.5 billions 
(dotted curve) and 2 billions (dashed curve), for this last value, the localization domain is 
located in the range [0.45 //, 0.72 //], it is larger compared to the case iV=1.5 billion ([0.45 //, 
0.60 //]). In figure |1~6|, we have fixed N at 1 billion and computed localization curves for R = 
10 fi (plain curve), i?=12.5 /a (dotted curve), and R — 15 \\i (dashed curve). In the last case, 
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the minimum of the localization curve is 1.1 for A = 5/i, but there is no more localization. 

The same effects occur in the case of a layered medium. When we increase the number 
of disks (figure |I7|) the localization domain enlarges, oppositely, when the thickness of the 



medium increases (figure [18]), the localization range narrows and disappears. 



B. A Numerical check 

In order to numerically control the behaviour of the electromagnetic field near a localiza- 
tion range, we have used a multidipolar diffusion formulation |L0| . The principle consists to 
discretize a dielectric sphere of radius a < A by an array of dipoles | IT[| , and to compute the 



scattered field using a multidipolar expansion. To characterize the scattered field, we use 



the intensity functions described in []12 |, where i and j respectively refer to the polar- 



ization of the incident and scattered electric field, 9 is the scattering angle. More precisely, 
we are interested by the intensity behaviour around the backscattering direction, where the 
localization phenomenon manifests by the creation of a backscattering peak. 

We compute the averaged intensity functions by supposing the disordered medium er- 
godic, which means that the mathematical expectation value of the intensity function is 
obtained by averaging intensity functions on a large number of disordered medium configu- 
rations. 

Due to the large number of spheres we used for computations, we simplify the problem 
by modeling a little sphere by a single dipole located at its center. 

C. Case of a spherical volume 

We first use the localization parameters previously obtained for the set of spheres in 
a spherical volume to compute the averaged intensity functions. In this model, a number 
of spheres around 1 million is too large to be handled by numerical computations. But 
we notice that the fundamental parameter which occurs in the localization criterion is the 
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N 

density n — — . To simulate an equivalent configuration, we define a medium of the same 
density but with a smaller number of spheres. However, the consequence is also to reduce 
the radius of the medium, and then alterate the macroscopic homogenization condition (34). 

We choose the set of the following parameters: R = 0.2/z and N = 1000. The localization 
parameters are then 

N = 1000 
a = 0.01 /Lt 

e s =!6 
R = 0.2 n 

X E [0.0805^,0.083^] 

We propose to compute the intensity functions for several values of N, when we approach 



a localization range (i.e., iV — > 1000). Figure [19] shows the intensity functions I xx (0) (in 
arbitray units) for the respective values of N =300,500,700,900. We notice the appearance 
of a peak in the backscattering direction when the value of N approaches the theoretical 
value of N given by the mathematical model (iV=1000). 



D. Case of a layered medium 

We use the localization parameters given for the set of spheres in a layered volume to 
compute the averaged intensity functions. Like in the previous example, the number of 
spheres is too large for a numerical computation. To overcome this difficulty, we use an 
equivalent medium with the same density. 

We choose the set of equivalent parameters: e x = 0.3//, e y = 0.3/x, e z = 0.28/i and 
iV = 800 . The localization parameters are then 

N = 800 
a = 0.01/i 
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e x = e y = 0.3/x 

e z = 0.28/1 

A E [0.081/4,0.083//] 

We present the intensity function curves I xx (6) when the number of spheres N have 
respectively the values 200, 400, 600, 800 (figure ^). We also observe a backscattering peak 
when N tends to the initial value where one observes an electromagnetic energy localization 
in the mathematical model. 



V. CONCLUSION AND PERSPECTIVES 

We have analytically calculated the value of the k e .l product as a function of the incident 

wavelength in order to find wavelength domains satisfying the localization criterion k e .l < 
r 3~ 

in different media. We have considered some configurations where the volume may be 

7T 

finite or infinite, and the scatterers are dielectric spheres or perfectly conducting little disks. 
We have studied the influence of electromagnetic system parameters on the localization 
domain, for instance, the complex permittivity of spheres, their radius, their number and 
the dimensions of the surrounding medium. Computations have shown that the density 
offers the optimal way to adjust the width of the localization domain. 

The multidipolar model allowed us to detect a backscattering peak when the system 
parameters are closed to the localization parameters provided with the theoretical model. 
Many works were already realized in order to calculate the expression of the backscattering 



peak [ffl, [13 1, however these results were not related to the localization criterion. This 
peak, predicted by the theoretical model, occurs from the crossed diagrams contribution 
responsible of the backscattering, it represents a first manifestation of the localization 
phenomenon. However, we were limited by the size of the matrix describing the system, and 
the conditions where the medium is macroscopically homogeneous were not exactly fulfilled. 
The same matrix limitation implies that we have, in the simulation, described a little sphere 
by a single dipole. A way to perform a more realistic simulation would consist to describe 
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the sphere itself by a set of dipoles located on a cubic lattice inside the sphere 



11 
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FIG. 2. Scheme for N spheres uniformly distributed in a spherical volume or in a layered volume 
of finite thickness. 
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FIG. 3. Scheme for N disks uniformly distributed in a spherical volume or in a layered volume 
of finite thickness. 
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FIG. 4. Localization curve for a set N = 1 million spheres of radius a = 0.01^, and relative 
permittivity e = 16, uniformly distributed in a finite spherical volume of radius R = 2fj,. 
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FIG. 5. Localization curves for a set of N = 1 million spheres, of relative permittivity e s = 16, 
and respective radii a = 0.095 (plain curve), a = 0.01 (dotted) and a = 0.0105/U (dashed), uniformly 
distributed in a finite spherical volume of radius R = 2/j,. 
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FIG. 6. Localization curves for a set of N = 1 million spheres, of radius a = 0.01/x, and respective 
relative permittivities e s = 15 (plain curve), e s = 16 (dotted) and e s = 17 (dashed), uniformly 
distributed in a finite spherical volume of radius R = 2/x. 



22 



kl 




0.075 0.080 0.085 



wavelength (micron) 

FIG. 7. Localization curves for a set of N = 1 million spheres, of radius a = 0.01/U, with a 
relative real permittivity e s = 16, and respective relative complex permittivities: (plain curve), 
0.4 (dotted), and 0.8 (dashed) and 1.2 (long dashed), uniformly distributed in a finite spherical 
volume of radius R = 2jx. 
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FIG. 8. Localization curves for the respectives sets of N = 1 million (plain curve), TV" = 1.5 mil- 
lions (dotted) and N = 2 millions (dashed) spheres, of radius a = 0.01/x, and relative permittivity 
e s = 16, uniformly distributed in a finite spherical volume of radius R = 2fi. 
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FIG. 9. Localization curves of a set of N = 1 million spheres of radius a = 0.01//, and relative 
permittivity e = 16, uniformly distributed in finite volumes of respective radii R = 2.1 (plain 
curve), 2.2 (dotted), 2.3 (dash) and 2.4 \x (long dashed). 
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FIG. 10. Localization curve of a set of N = 10 5 spheres of relative permittivity e s = 16, and 
radius a = O.Ol/i uniformly distributed in a layered volume of thickness e z = 1.5^, and transverse 
lengths e x = 1.5// and e y = 1.5/t. 
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FIG. 11. Localization curves for the set of N = 10 5 (plain curve), TV = 1.210 5 (dotted) and 1.510 5 
(dashed) spheres of radius a = 0.01//, and relative permittivity e s = 16, uniformly distributed in a 
medium of thickness e z = 1.5// and of transverse lengths e x = 1.5/t and e y = 1.5//. 
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FIG. 12. Localization curves for a set of iV = 10 5 spheres of radius a = 0.01// and relative 
permittivity e s = 16, uniformly distributed in different media of respective thickness e z = 1// 
(plain curve), e z = 1.5// (dotted) and e z = 1.8// (dashed), and of respective lengths e x = 1.5// and 
y = 1.5/t. 
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FIG. 13. Localization curve for a set of N = 1 billion of perfectly conducting disks of radius 
= 0.05//, uniformly distributed in spherical volume of radius R = 12.5/t. 



kl 




wavelength (micron) 



FIG. 14. Localization curve for a set of iV=150 millions perfectly conducting disks, of radius 
= 0.05// uniformly distributed in a layered volume of thickness e z = 10// and transverse widths 
; = 10// and e y = 10//. 
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FIG. 15. Localization curves for the respective sets of N=l (plain curve), N=l.5 (dotted), N=2 
billions (dashed) perfectly conducting disks of radius a = 0.05/U, uniformly distributed in a spherical 
volume of radius R = I5fi. 
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FIG. 16. Localization curves for a set of N=l billion disks of radius a = 0.05/i uniformly 
distributed in different media of respective radii R = 10 (plain curve), i?=12.5 (dotted), R = 15fj, 
(dashed). 
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FIG. 17. Localization curves for the respective sets of iV=150 (plain curve), 7V=200 (dotted) 
and iV=250 millions (dashed) perfectly conducting disks of radius a = 0.05//, uniformly distributed 
in a layered medium of thickness e z = 10/i, and transverse widths e x = 10// and e y = 10//. 
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FIG. 18. Localization curves for a set of N = 150 millions disks of radius a = 0.05//, uniformly 
distributed in several media of respective thickness e z = 10// (plain curve), e z = 15// (dotted) and 
e z = 20// (dashed), and transverse widths e x = 10// and e y = 10//. 
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FIG. 19. Averaged intensities over 20 configurations as a function of the scattering angle 9, for 
an x incident polarization of the electric field for different media, where N =300, 500, 700, 900 
spherical scatterers of radius a = 0.01^, and relative permittivity e s = 16 uniformly distributed in 
a spherical volume of radius R = 2/x. The incident wavelength is Aq = 0.082^. 
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FIG. 20. Averaged intensities over 20 configurations as a function of the scattering angle 8 for 
an x incident polarization of the electric field in a layered medium of thickness z = 0.28/x, where 
iV=200, 400, 600, 800 spherical scatterers of radius a = 0.01/x, and relative permittivity e s = 16. 
The incident wavelength is Aq = 0.082^. 
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